AD-A268  929 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


I 


r’  iDi'c  'PDOrTK'.j  Doraen  for  thi^  coHection  of  .rifofrraiton  <v  estimated  ro  average  ■  pouf  oer  '«por\e.  ncluding  the  time  *or  reviewing  mttructionv,  searching  evistmg  data  tourcn.  1 
jatherirtg  ind  mamtainmg  the  data  needed,  and  camoteting  and  rev«ew*rg  the  collection  of  information  Send  comments  regarding  this  Ouroen  estimate  or  any  other  asoect  of  this 
cciiection  0*  'ntormaiion.  tncsuding  suggestions  tor  reducing  this  ourden,  to  Washington  Headguarters  Services.  Oireaorate  tor  information  Operations  and  Reports.  ’2 15  Jefferson 
Oawis  M  gnwav,  suite  1204.  Arlington.  VA  22202-4302.  and  to  the  Office  of  Manageme"!  and  Budget.  Paperwork  Reduaion  Project  (0704-0198),  Washington.  DC  20503 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

1/11/93 


3.  REPORT  TYPE  ANO  OATES  COVERED 

FINAL  15  Nov.  91  -  lA  Jfov. 


92 


4.  TITLE  AND  SUBTITLE 

Fundamental  Research  in 
Aerospace  Applications 


Nonlinear  Random  Vibration  for 


6.  AUTHOR(S) 


Y.  K.  Lin 


AFDSfi-TS- 


5.  FUNDING  NUMBERS 


G-AFOSR-91-0073 


0  6  6  6 


7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  AODRESS(ES) 

Center  for  Applied  Stochastics  Research 

College  of  Engineering 

Florida  Atlantic  University 

500  N.W.  20th  Street 

Boca  Raton.  FL  33431-0991 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  ANO  AOORESS(ES) 


AFOSR 

Building  410 
Bolling  AFB 

Washington,  0.  C.  20332^^448^ 


10.  Sff)NjOR|lG/t4flNITORING 
AeMCY  KPCfT%UMBER 


ELECTE 
AUG  3  1 1993 


11.  SUPPLEMENTARY  NOTES 


A 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 


Ajp.'ovod  tor  r-itl  lo  rolenaail 

O  dlatPltutlonuiilieltBA 

8  30  Ong 

13.  ABSTRACT  (Maximum  200  words)  ^ 


q3-20243 

_ inn  uni  HI  111 


Response  and  safety  of  nonlinear  structures  under  random  excitations  are  investigated.  The  excitations 
can  be  either  external,  or  parametric,  or  both.  Extensions  are  made  along  two  directions  from  the 
earlier  works  on  exact  solutions  when  the  excitations  are  Gaussian  white  noises:  (1)  the  excitations 
are  impulsive,  but  non-Gaussian,  (2)  the  system  does  not  belong  to  the  class  of  generalized  stationary 
potential,  thus,  exact  solution  is  not  obtainable  at  the  present  time.  In  the  first  case,  the  probability 
density  of  the  structural  response  is  governed  by  a  partial  differential  equation,  with  infinite  number 
of  terms.  A  perturbation  analysis  is  devised  to  obtain  approximate  solutions.  In  the  second  case,  a 
new  technique  is  developed  to  obtain  approximate  solutions  on  the  principle  of  weighted  residuals. 
On  the  issues  of  structural  safety,  two  failure  modes  are  investigated:  catastrophic  failure  and  fatigue 
failure.  Catastrophic  failure  occurs  when  structural  response  exceeds,  for  the  first  time,  a  critical 
limit.  Fatigue  failure  is  investigated  from  a  fracture  mechanics  pointAjf  view,  namely,  fiulure  is 
assumed  to  occur  when  a  dominant  crack  is  propagated  to  a  critical  size. 


14.  SUBJEQT  TERMS 

Random  Vibration,  Nonlinear  Systems 

15.  NUMBER  OF  PAGES 

21 

16.  PRICE  CODE 

1  17.  SECURITY  CLASSIFICATION 

1  OF  REPORT 

j  Unclassified 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

Unclassified 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

Unclassified 

20.  LIMITATION  OF  ABSTRACT 

UL 

INTRODUCTION 

Aerospace  structures  are  often  subjected  to  random  excitations.  One  major  source  of 
random  excitations  is  the  air  turbulence,  which  may  be  generated  by  natural  geo-thermal 
processes,  afflux  of  a  propulsion  system,  flight  maneuver,  or  explosions  in  a  battle  field. 
Depending  on  the  way  in  which  such  random  excitations  interact  with  a  structure,  they  can  be 
classified  into  two  types:  the  external  (or  additive)  excitations  and  the  parametric  (or 
multiplicative)  excitations.  In  the  equations  of  motion,  external  excitations  appear  as 
inhomogeneous  terms,  whereas  parametric  excitations  appear  in  the  coefficients  of  the 
unknowns. 

Aerospace  structures  are  often  operated  into  the  nonlinear  regime  due  to  elevated 
temperature  or  large  deformation.  Furthermore,  certain  fluid-structure  interaction  phenomena 
are  known  to  be  highly  nonlinear,  e.g.,  stall  flutter,  vortex  induced  vibration,  and  high  angle 
of  attack  during  maneuver. 

The  objective  of  the  research  project  is  to  develop  solution  techniques  for  nonlinear 
aerospace  structures  under  either  additive  or  multiplicative  random  excitations,  or  both.  The 
required  solutions  include  probability  densities  or  statistical  properties  of  structural  response  to 
random  excitations,  as  well  as  reliability  of  such  structures  to  fulfill  their  intended  functions. 


=  ffx,  /) 


j  =  1,  2,  N\ 


it  =  1,  2,  ...,  M  (1) 


where  fj  and  gj^  are  nonlinear  functions,  and  W^Ct)  are  random  excitations.  The  possibility  for 
exact  solutions  exists  only  when  Wk(t)  are  Gaussian  white  noises.  In  this  case  X(t)  is  a  Markov 
vector,  whose  probability  density  p(x)  at  the  state  of  statistical  stationarity  is  governed  by  the 
following  reduced  Fokker-Planck  equation  [e.g.l] 


where 


(2) 


G,  =  4(x)p(jr)  -  I  ^[Bj,(x)pix)] 


(3) 


Ajix)  =fj(x)  +  gj^(x)  (4) 

8jt  8ks 

and  where  K,,  is  the  cross-spectral  density  of  W,(t)  and  W,(t).  In  Eqs.  (2)  -  (5),  a  lower  case 
symbol  represents  a  possible  value  of  the  corresponding  upper  case  random  quantity.  The 
second  order  partial  differential  equation,  Eq.  (2),  is  solved  with  the  condition  Gj  =  0  at  all  the 
boundaries.  We  have  obtained  previously  solutions  for  a  class  of  nonlinear  systems,  called  the 
class  of  generalized  stationary  potential  [2-4]. 

The  present  effort  is  directed  at  relaxing  some  idealized  assumptions.  In  some 
aeronautical  applications,  the  excitations  on  a  structure  are  sometimes  irregular  intermittent  short 
pulses,  which  can  be  reasonably  modeled  as  a  white  noise;  however,  it  may  be  non-Gaussian. 
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In  this  case,  Eq.(2)  is  replaced  by  the  following  partial  differential  equation  of  an  infinite  order. 


2!  dx.dx, 


3!  dx,dXjdx^ 


1  3^ 

4!  ^x^^XJ^x^^x, 


=  0 


(6) 


where  A;  and  Bjj  have  been  defined  in  Eqs.  (3)  and  (4),  whereas  Cijk,  Djju,  etc.,  are  calculated 
according  to 


lim  —  £[AX.AX,AXjAr(r)  =x] 

4/-.0  At  ^ 


(7) 


D.  =  lim  —  £[AXAX;.AX*AX,|;!f(0  =x] 

■'  A/-.0  At 


(8) 


In  Eqs.  (7)  and  (8),  AXj  =  X^t+At)  -  Xi(t),  etc.  A  closed  form  solution  for  Eq.  (6)  is  clearly 
impossible. 

A  perturbation  scheme  is,  therefore,  devised  to  obtain  an  approximate  solution  for  non- 
Gaussian  white  noise  excitations.  For  simplicity,  consider  the  case  of  a  single  white  noise 
excitation,  which  may  be  treated  as  a  sequence  of  independently  arriving  impulses,  each  with 
a  random  magnitude  Y.  The  probability  density  of  Y  need  not  be  symmetrical.  If  the  average 
arrival  rate  of  the  random  impulses  is  \,  then  the  departure  from  Gaussianity  can  be  measured 
by  the  products  X  EfY""^^]  where  n  are  integers.  A  Gaussian  white  noise  rqpresents  the  limiting 
case  X  E[Y“^^]  -►  0  for  any  n  >  0.  To  account  for  non-Gaussianity,  let 
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X£[y»-2]  =  e-  n  =  0,  1,  2,  ...  (9) 

Then  an  approximate  solution  to  Eq.(6)  may  be  taken  in  the  form  of  [S] 

P(x)  =  Po(Jc)  [1  +  eQ^(x)  +  e^QjCx)  +  ...]  (10) 

Detailed  solutions  for  Q,(x)  and  Q2  are  given  in  [5]. 

Two  specific  nonlinear  systems,  under  external  and  parametric  excitations  respectively, 
are  investigated  in  detail.  Typical  results  are  illustrated  in  Figs.  1  and  2,  These  analytical 
results  compared  favorably  with  the  Monte  Carlo  simulation  results.  It  can  be  seen  that  the 
effect  of  non-Gaussianity  in  the  excitation  is  significant  for  the  externally  excited  system,  but 
negligible  for  the  parametrically  excited  system,  even  though  the  two  systems  share  the  same 
probability  distribution  when  the  impulsive  excitation  is  Gaussian.  Moreover,  increasing  level 
of  nonlinearity  in  the  parametrically  excited  system  tends  to  reduce  the  effect  of  excitation  non- 
Gaussianity,  whereas  the  opposite  is  true  for  the  externally  excited  system.  For  practical 
purposes,  an  external  impulsive  noise  excitation  may  be  treated  as  a  Gaussian  white  noise  if  the 
product  of  the  average  pulse  arrival  rate  and  the  relaxation  time  of  the  dynamical  system  is  about 
20  or  larger  [5]. 

METHOD  OF  WEIGHTED  RESIDUAL 

When  a  system  does  not  belong  to  the  class  of  generalized  stationary  potential,  an  exact 
probabilistic  solution  is  presently  unobtainable.  Various  approximation  techniques  have  been 
developed  in  the  past,  the  most  recent  one  being  the  dissipation  energy  balancing  technique 
proposed  by  Cai  and  Lin  [6].  To  explain  this  particular  technique,  consider  a  second  order 
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nonlinear  system,  governed  by 


The  reduced  Fokker-Planck  equation  associated  with  system  (10)  is  given  by 


(11) 


X2 -^p  *  (JC,  ,x^)  ^  [ -Hix^  ,x^)  *  xXy^(x,  .Xj)  ]p  •  (Xj  } 

OJC^  0  Xf^ 

-xX  J1  [/(X,  ,X2)/(X,  ,X^)p  *  (X,  .Xj)]  =  0 
3x2 


(12) 


where  p*  (x^,x^)  is  the  stationary  probability  density  of  the  response  for  system  (11),  Let  us 


suppose  that  the  oscillator  does  not  belong  to  the  class  of  generalized  stationary  potential, 
therefore,  an  exact  solution  p  *  (x,  ,x^)  for  Eq.(12)  is  presently  not  obtainable.  Since  the  class 

of  generalized  stationary  potential  is  the  broadest  solvable  class  known  to-date  which  includes 
a  large  class  of  linear  and  nonlinear  systems,  a  logical  approximation  scheme  is  to  replace  the 
original  system  (11)  by  an  equivalent  nonlinear  system  within  the  class  of  generalized  stationary 
potential,  say 


X  *  A(X,X)  (13) 

Note  that  the  right  hand  sides  of  Eqs.(ll)  and  (13)  can  be  made  the  same.  The  strategy  is  to 
find  an  h(X,X)  function  within  the  class  of  generalized  stationary  potential  which  is  closest  to 

the  original  H{X,X)  function.  The  criterion  proposed  in  [6]  for  selecting  such  an  h{X,^ 

function  is  that  the  average  dissipation  energy  remains  the  same  for  the  substituting  and  the 

substituted  systems. 

« 

The  dissipation  energy  balancing  criterion  has  been  generalized  to  the  criteria  of  zero 
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weighted  residual  in  [7],  It  can  be  shown  [7]  that  such  criteria  are  equivalent  to 

=£{[A(X,,X,)-H(X,,Xi))£;^^^)  =0  (14) 

where  X,  =  X,  JXj  =  X  and  M(Xp  X^)  is  a  suitable  weighing  function  of  X,  and  Xj,  and 

£[  ]  denotes  the  ensemble  averaging  with  respect  to  the  approximate  probability  density 
p(Xj  ,^2) .  The  residual  depends  on  the  choice  of  M.  By  judiciously  selecting  the  weighing 
functions,  the  approximation  can  be  made  closer  and  closer  to  the  exact  probability  density.  The 
choice  of  M  =  X2  is  equivalent  to  the  requirement  of  dissipation  energy  balancing,  which 
remains  the  most  appropriate  and  crucial.  The  remaining  criteria  proposed  in  [7]  are  of  the  form 
Af  =  XjX*  where  k  =  0,  1,  2,  ...,  n. 

It  is  of  interest  to  note  that  the  well-known  procedure  of  equivalent  linearization  [8]  is 
a  special  case  of  the  present  procedure,  in  which  h(Xi,  is  assumed  to  be  linear,  namely 
h(Xi,X2)  =  aX,  +  /SXj.  However,  the  equivalent  linearization  procedure  is  unsuitable  to  treat 
parametric  excitations;  thus,  its  applicability  is  restricted  to  the  case  of  external  excitations, 
namely,  when  f^  functions  in  JEq.(U)  are  constants.  Another  special  case  is  the  procedure  of 
equivalent  nonlinear  systems  [9],  proposed  also  for  the  case  of  external  random  excitations,  but 
allowing  the  H  function  in  Eq.  (11)  to  be  nonlinear.  It  is  expected  that  neither  of  the  special 
cases  is  as  accurate  as  the  method  of  weighted  residual,  because  in  each  of  these  procedures  the 
approximation  is  selected  only  from  a  sub-class  within  the  much  larger  exactly  solvable  class  of 
generalized  stationary  potential. 

Approximate  probability  densities  obtained  from  three  different  procedures,  namely. 
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equivalent  linearization,  equivalent  nonlinear  system,  and  weighted  residual,  are  compared  in 
Figs.  3  and  4.  It  is  seen  that  the  method  of  weighted  residuals  is,  indeed,  consistently  superior. 


COMBINED  HARMONIC  AND  RANDOM  EXCITATIONS 
Consider  the  following  equation  of  motion 

X  +  2^0}  *  JoX  *  hiX,X)  =^.(X,X)sin(/V0  +  1(0 

where  |(t)  is  a  broad-band  random  process,  playing  the  role  of  external  excitation.  The  system 

is  also  excited  by  sinusoidal  parametric  excitations  sinOj/t);  j  =  1,2 . N.  A  sinusoidal 

excitation  sin(/V0  is  important  only  when  near  tuning  occurs;  namely,  when  |  jv  -  \  is 

small.  Otherwise,  the  presence  of  the  sinusoidal  parametric  excitations  can  be  ignored,  and  the 
problem  is  reduced  to  one  with  the  broad-band  external  excitation  alone. 

In  the  near  tuning  case,  the  stochastic  averaging  technique  fails  to  de-couple  the 
amplitude  process  from  the  phase  process.  However,  an  exact  solution  procedure  is  devised  for 
a  class  of  nonlinear  systems.  When  exact  solutions  are  not  obtainable,  approximate  solutions 
can  be  obtained  by  applying  the  method  of  weighted  residuals.  The  accuracy  of  the  approximate 
results  is  substantiated  by  Monte  Carlo  simulations.  A  technical  paper  has  been  submitted  for 
journal  publication  [10].  Typical  results  are  shown  in  Fig.  5. 

CATASTROPHIC  FAILURE 

Catastrophic  failure  occurs  when  a  response  quantity  of  interest  exceeds,  for  the  first 

time,  a  critical  magnitude.  The  time  at  which  a  catastrophic  failure  occurs  is  known  as  the  first- 
« 

passage  time  in  the  theory  of  stochastic  process.  For  aerospace  structures,  the  most  relevant 
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quantity  for  the  catastrophic  failure  analysis  can  be  either  the  response  amplitude,  or  total  energy 
(the  sum  of  potential  and  kinetic  energies).  Under  certain  conditions,  such  a  critical  quantity 
may  be  modeled  as  a  Markov  process.  The  critical  state  is  then  represented  by  the  absorbing 
boundary  of  the  process.  The  statistical  moment  =  EU"]  of  the  failure  time  T  is  governed 
by  the  recursive  equation 

ia2(Zo)  /x„ . ,  +  n  =  0,1,2,...  (16) 

2  dzS  ^^0 

where  Zq  is  the  initial  state,  which  is  treated  as  the  independent  variable.  By  definition,  /io  = 
1.  The  coefficients  m(zo)  and  o^(zo)  are  called  the  drift  and  diffusion  coefficients,  and  are 
determined  from  the  equations  of  motion,  and  the  spectral  densities  of  the  random  excitations. 
Equation  (16)  is  solved  subject  to  the  boundary  conditions 


l^n  *  1  (^c)  =  0 

M„*i(0)  =finne 


(17) 


where  is  the  failure  state.  In  the  special  case  n  =  1,  Eq.(17)  is  reduced  to  the  well-known 
Pontryagin  equation  [11].  It  is  shown  that,  when  diffusion  coefficient  vanishes  at  a  boundary, 
known  as  a  singular  boundary,  it  must  be  either  an  entrance  or  a  regular  boundary  in  the  sense 
of  Fuller,  in  order  that  the  solution  is  physically  meaningful. 

A  general  method  of  solution  for  Eq.(16)  has  been  devised  for  the  following  nonlinear 
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system 


^  *  e/i(X,AO  >  (4X  =  e^gjiX,X)^.it) 

where  |j(t)  are  broad-band  noises.  A  technical  paper  documenting  the  detailed  analysis  has  been 
accepted  for  publication  in  the  Journal  of  Applied  Mechanics  [12].  Computed  results  for  a 
Duffing  oscillator  is  shown  in  Fig.  6. 

FATIGUE  CRACK  GROWTH 

From  a  fracture  mechanics  point  of  view,  fatigue  failure  may  be  treated  as  propagation 
of  a  dominant  crack  to  a  critical  size.  Our  study  is  based  on  the  following  crack  growth  model, 

(19) 

at 

where  A(t)  is  the  size  of  a  dominant  crack,  g(  )  represents  a  crack  propagation  law,  for 
example,  the  well-known  Paris-Erdogan  law,  Y(t)  accounts  for  the  random  variation  of  material 
resistance,  fi  is  the  average  number  of  stress  maxima  per  unit  time,  and  AX  is  the  random  stress 
range  defined  as 

AX(r)  =  |X(r,)  1  ,  t,  <  t  <  t,  (20) 

and  where  t,  and  t2  are  the  times  at  which  two  neighboring  extrema  of  random  stress  X(t)  occur. 
A  procedure  is  developed  to  calculate  the  probability  of  fatigue  failure,  namely,  for  the  dominant 
crack  to  grow  to  a  critical  crack  size  [13).  This  procedure  is  applicable  when  the  fatigue  crack 

t 

growth  is  a  slow  process  compared  with  the  random  stress  process,  which  is  the  case  with  the 
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high  cycle  fatigue.  In  the  special  case  in  which  the  stress  is  a  narrow-band  stationary  Gaussian 
process,  analytical  expressions  have  been  obtained  for  the  probability  densities  of  the  crack  size, 
the  fatigue  life,  and  the  liability  function. 

Calculations  have  been  carried  out  for  the  case  of  a  central  crack,  depicted  in  Fig.  7. 
The  computed  probability  density  is  illustrated  in  Fig.  8,  along  with  Monte-Carlo  simulation 
results. 

CONCLUDING  REMARKS 

In  this  report,  we  have  high-lighted  the  research  accomplishments  during  the  two  years 
under  grant  G-AFOSR-9 1-0073.  A  total  of  ten  publications  have  been  generated,  with 
acknowledgement  to  the  AFOSR  support.  It  is  hoped  that  the  information  provided  in  these 
highlights  will  direct  an  interested  reader  to  more  detailed  accounts  in  the  quoted  references. 
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1  •  Stationary  probability  density  of  displacement  X  for  a  Duffing  oscillator  under 

external  impulsive  excitations.  Excitation  pulse  amplitude  non-symmetrically 
distributed. 
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Stationary  probability  density  of  displacement  X  for  a  nonlinear  system  under 
parametric  impulsive  excitations.  Excitation  pulse  amplitude  non-symmetrically 


distributed. 


Fig.  3 


Fig.  4. 
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Mean  square  displacement  for  system  X  +  +  yX^  +  5X^  ~  W{t) 

0  =  0.1,  5  =  0.5 


Mean  square  displacement  for  system  X  0X  +  yX^  +  6X^  =  W(t) 


jS  =  0.01,  5  =  1.5 


Stationary  probability  density  of  amplitude  A  for  nonlinear  system 

X  +  2i:uoX  *  +  (4X[l  *  Xsin(2;/0]  =  W(t), 


j8  =  0.15,  r  =  0.05,  Wo  =  10,  V  =  11,  X  =  0.1 
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